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We present high statistics simulations of weighted lattice bond animals and lattice trees on the 
square lattice, with fugacities for each non-bonded contact and for each bond between two neigh- 
bouring monomers. The simulations are performed using a newly developed sequential sampling 
method with resampling, very similar to the pruned-enriched Rosenbluth method (PERM) used 
for linear chain polymers. We determine with high precision the line of second order transitions 
from an extended to a collapsed phase in the resulting 2-dimensional phase diagram. This line 
includes critical bond percolation as a multicritical point, and we verify that this point divides the 
line into different universality classes. One of them corresponds to the collapse driven by contacts 
and includes the collapse of (weakly embeddable) trees. There is some evidence that the other is 
subdivided again into two parts with different universality classes. One of these (at the far side from 
collapsing trees) is bond driven and is represented by the Derrida-Herrmann model of animals having 
bonds only (no contacts). Between the critical percolation point and this bond driven collapse seems 
to be an intermediate regime, whose other end point is a multicritical point P* where a transition 
line between two collapsed phases (one bond-driven and the other contact-driven) sparks off. This 
point P* seems to be attractive (in RG sense) from the side of the intermediate regime, so there 
are four universality classes on the transition line (collapsing trees, critical percolation, intermedi- 
ate regime, and Derrida-Herrmann). We obtain very precise estimates for all critical exponents for 
collapsing trees. It is already harder to estimate the critical exponents for the intermediate regime. 
Finally, it is very difficult to obtain with our method good estimates of the critical parameters of 
the Derrida-Herrmann universality class. As regards the bond-driven to contact-driven transition 
in the collapsed phase, we have some evidence for its existence and rough location, but no precise 
estimates of critical exponents. 



I. INTRODUCTION 

Lattice animals are just clusters of connected sites 
on a regular lattice. They play an important role in 
many models of statistical physics, as e.g. percolation Q 
and the Ising model, where they are related to Fortuin- 
Kastleyn clusters and to the famous Swendsen-Wang al- 
gorithm [2J.LJ. Apart from this, they form also the pro- 
totype model for randomly branched polymers, just as 
self avoiding random walks (SAWs) are a model for un- 
branched polymers Q. 

In the case of single unbranched polymers in a very 
diluted solvent, a much studied phase transition happens 
when the solvent deteriorates (as it usually does when 
temperature is lowered). Below the so-called #-point a 
polymer no longer forms a swollen coil with Flory radius 
Rn ~ N v with v > 1/2 (N is here and in the follow- 
ing the number of monomers), but rather a collapsed 
"globule" with Rn ~ N 1 ^, where d is the dimensional- 
ity of space |E 01- ln the simplest model of polymers 
living on a regular lattice, the collapse is induced by 
an effective attractive interaction between non-bonded 
monomers on neighbouring lattice sites. It is easy to see 
that such a monomer-monomer attraction is equivalent 
to a monomer-solvent repulsion, so there is no need to 
include the latter, if one is only interested in universal 
properties of the transition. 

A similar collapse transition is expected to occur also 
for branched polymers. Here the situation is somewhat 
more complicated, though. In a lattice model, one can 
introduce two different kinds of attraction, so one obtains 



a two dimensional phase diagram with a line of collapse 
transition points. The basic reason for this difference 
to unbranched polymers (interacting self avoiding walks, 
ISAW) is that the number of bonds is fixed to b = N — 1 
for SAWs (as it is also for trees), while it can fluctuate, 
b > N — 1, for general animals. Thus one can introduce 
two different fugacities for monomer- monomer bonds and 
for non-bonded monomer-monomer contacts. As for un- 
branched polymers, there is no need for introducing a 
separate monomer-solvent interaction, since the number 
s of monomer-solvent contacts is not independent, but is 
given by 

AfN = 26 + 2fc + s, (1) 

where N is the lattice coordination number (N = 2d on 
a simple hypercubic lattice in d dimensions) and k is the 
number of non-bonded monomer- monomer contacts. A 
general partition sum for interacting lattice animals is 
therefore 0, H H E3, IH E3, EE E3, Il3 

z N (y,r) = J2 c Nbky h - N+1 r k , (2) 

b,k 

where Cjvbfc is just the number of configurations (up to 
translations and rotations) of connected clusters with N 
sites, b bonds, and k contacts. Notice that we changed 
the definition slightly with respect to 0, B H E [H 111 
If3.IT3.ITH. so that Zn(u, t) is non-zero for y = 0. 

The phase diagram in terms of the control parameters 
y and r is shown in Fig. 1 for 2-d animals on the square 
lattice. There is an extended phase for small y and r 
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FIG. 1: Phase diagram for interacting animals. The full line 
separates an extended phase (below) from a collapsed phase 
(above). The points on this line are obtained by Monte Carlo 
simulations in the present work, except for the point at y = 
6.48, r = which is taken from 16]. At y — the clusters are 
trees (minimal number of bonds) , while at r = they have no 
contacts but only bonds. The dashed line corresponds to bond 
percolation, with the critical point being at y — r = 2. The 
short dashed-dotted line is a rough estimate for the transition 
between a contact-rich and a bond-rich collapsed phase. 



which includes also the unweighted animal model (y = 
t = 1), and at least one collapsed phase. At the collapse 
transition we expect that 

Z N (y,T = T c (y))~rty) N N- 9 (3) 

where n(y) should depend continuously on y, but 8 
should take discrete values depending on the respective 
universality class. 

There have been claims, based on exact enumerations 
of very small clusters jjj , that there are two distinct 
collapsed phases, one bond-rich and the other contact- 
rich. This has led to some controversy, since later authors 
could not find it either with other numerical models |8( 
or in simplified models Q. We also find such a tran- 
sition between two collapsed phases, although at signif- 
icantly smaller values of y (short dashed-dotted curve). 
The two end points of the transition curve shown in Fig. 1 
are given by collapsing weakly embeddable trees (y = 0) 
and by a model where non-bonded contacts are forbidden 
(t = 0). The latter has been studied in detail by Derrida 
and Herrmann , and others pjj Oil HE • 

Bond percolation is a special model of weighted clus- 
ters. The generating function for bond percolation clus- 
ters is defined as 

z p N eic (p) = J2 c ^p b ( 1 -p^ k+s - ( 4 ) 

b,k 

More precisely, this equation gives the probability that 
the origin is connected to a cluster of precisely N sites, 
if lattice bonds are established with probability p and 



broken with probability 1 — p. Using Eq. QJ, one sees 
that 

Z»™(p) = (l-p)« N -™+ 2 p N -iZ N (y(p),T(p)), (5) 
with the curve (y(p),T(p)) parametrized as 

y=p/(l-p) 2 , T = l/(l-p), 0<p<l. (6) 
or given explicitly by 

y = t(t - 1) . (7) 

This curve is also shown in Fig. 1. It contains in par- 
ticular the critical bond percolation point, p — 1/2, cor- 
responding to y — t = 2. This point lies also on the 
collapse transition line Q, and divides it into two parts: 
One branch with high density of contacts (y < 2) and one 
with high density of bonds (y > 2). Since it is known 
that the critical percolation point is f ully repulsive in 
the renormalization group (RG) sense [21|. it would be 
natural to assume that the RG flow along the transition 
curve goes to the two end points y = (collapsing trees) 
resp. t — (Derrida-Herrmann model). Thus we could 
naively expect that the critical behaviour on each of the 
two parts of the transition curve is given by the fixed 
point at its end. But if there are indeed two different col- 
lapsed phases, then the point at which the transition line 
between them meets the collapse transition line should 
be some multicritical point (y*,r*). A priori, this point 
could be the critical percolation point. But according 
to 0, El this is clearly not the case: while y = 2 
for percolation, y* ss 4 (see Fig. 2 of pj). If all that is 
true, then there should be two multicritical points on the 
transition curve, and three different universality classes 
in addition to the percolation class. 

In the following we shall investigate the details of this 
scenario by means of extensive Monte Carlo simulations. 
In Sec. 2 we shall shortly describe the algorithm. The 
critical percolation point is studied in Sec. 3, both as a 
means to verify the exactness of our codes and to dis- 
cuss the cross-over exponents at this point. We shall 
see that both cross-over exponents (the one associated to 
fluctuations of the number of bonds and the other one 
for contacts) are exactly equal to 1/2, although the fluc- 
tuations of the number of contacts has huge corrections 
to scaling. After that we turn in Sec. 4 to the collapse 
of trees, i.e. to the case y < 2. We find that there 
our results are very clean, allowing us to obtain conjec- 
tures for the exact values of all critical exponents. We 
also verify with high precision a previous conjecture that 
the transition curve for y < 2 is strictly horizontal, with 
r = 2. Next, we will turn to the case y > 2 in Sec. 5. 
There we have much more problems in arriving at a con- 
sistent scenario, presumably because there exist indeed 
two different collapsed phases which meet on the transi- 
tion curve at a point different from critical percolation. 
In particular, at face value our results would suggest that 
the cross-over exponent for bonds depends continuously 
on y, and would thus be not universal. But apart from 
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these problems it seems very clear that the critical be- 
haviour for y > 2 is not in the same universality class as 
that for y < 2, in contrast to recent speculations [14j . We 
follow up the question of a transition between two col- 
lapsed phases by simulating deep in the collapsed phase. 
We conclude with a discussion in Sec. 7. 



small samples and will have rather large errors. With re- 
sampling we can work at p < p c , and we will have many 
more large finished clusters in the sample. 

All the work reported in this paper was done on fast 
Linux PCs and used nearly 2 years of CPU time. 



III. THE PERCOLATION POINT 



II. NUMERICAL METHODS 

We use essentially the same method as in |23) . 
There we introduced a sequential sampling method 
with re-sampling, implemented depth-first as in the 
PERM (pruned-enriched Rosenbluth method) algorithm 
of Ref . [2J| . More precisely, we first choose a value of p 
near the percolation threshold (i.e., p ps 1/2), and sim- 
ulate bond percolation clusters at this p by means of a 
variant of the Leath algorithm [2j|. While they are still 
growing we estimate their current contribution to the 
percolation partition sum, and - by re-weighting them 
according to the animal ensemble, Eq. (J^Jl - to the an- 
imal ensemble. If the latter happens to be higher that 
average, we clone the cluster and let both copies evolve 
independently further. If the weight is too small, we 
prune with probability 1/2 and increase the weights of 
the survivors by a factor 2. 

For this to be efficient, we have to take the following 
considerations into account: 

- The clusters are grown breadth first, not depth 
first |2^| , although the "population control" is done depth 
first. 

- The optimal value of p depends on y and r, and has to 
be found by trial and error. 

- The optimal value of the pruning/cloning threshold was 
found, to very goo d precision, to be the same as for or- 
dinary animals 23] , although we could improve also this 
by trial and error search. 

- To make things even more subtle, the optimal param- 
eters depended slightly on the maximal sizes to be simu- 
lated. 

To monitor the performance of the algorithm, we used 
mostly the variance of the estimated partition sum. No- 
tice that estimating the partition sum is an integral part 
of this algorithm (it is needed for the "population con- 
trol"). In addition, we checked that the distributions of 
"tour weights" remain acceptable, as described in |23J. 

It is clear that this algorithm should work best near the 
percolation point (since there is then a minimal amount 
of re-sampling needed). Indeed, exactly at the perco- 
lation point the method works even better than Leath 
itself, if one uses p slightly below 1/2. The reason is 
that with p = 1/2 and without re-sampling the sample 
contains very few large but finite clusters. The parti- 
tion sum Eq. (QJ decreases at p = p c = 1/2 very slowly 
with N, like A/- - 054945 pj. Thus most clusters which did 
not stop growing early will do so only very late, beyond 
the maximal N we can simulate. Therefore, estimates of 
specific heat or of gyration radii will be based on rather 



We use the above algorithm to simulate critical bond 
percolation clusters (y = t = 2 } p = 1/2) of up to 
N = 6000 sites. The nominal value at which the 
clusters were grown was p ~ 0.473, the difference be- 
tween this and the target value p — 1/2 being made up 
by re-weighting. In this way the variances of the observ- 
ables were reduced by roughly two orders of magnitude, 
as compared to a straightforward Leath algorithm. 

Results for the logarithm of the partition sum and for 
the gyration radius will be shown later, together with 
results obtained for y < 2. They will be discussed in 
detail in the next section, here we just point out that 
they are in perfect agreement with the expectations [l| 



Z p N eTC {p= 1/2) ~N-^ 91 



R 



N 



^96/91 (8) 



and with simulations using the plain Leath algorithm. 

More interesting are the average numbers of contacts 
and bonds, and the fluctuations thereof. Near the critical 
point, one has the scaling ansatz p| 



Z V n TC (p) ~ N- 5 / gl F{(p - l/2)N a ) 



(9) 



with a = 36/91. This gives <91n Z^ crc /dp\ p=1/2 ~ N a 
and <9 2 lnZP cr Wl P =i/2 ~ N 2a . On the other hand, 
using the defining equation (3), one obtains at p = 1/2 



d\nZ p N crc 



dp 



p=l/2 



= -(b-8-k) 



(10) 



and 



d 2 \nZ p N CIC 



dp 2 



p=l/2 



= ~{Vax[b-s-k]-(b+s + k)} . (11) 



Using this and eliminating s in favour of N, one finds 
that 



while 



(36 + fc) =4N + 0(N a ) 



Var[36 + /c] = 2(6) + 0(N 2a ) 



(12) 



(13) 



Our estimates for the average numbers of bonds and 
contacts are 



(b) = 1.1215(3)7V + 0(7V CT ), (k) = 0.6360(5)JV + O(7V° 



in excellent agreement with Eq. I|12fl, 



(14) 
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FIG. 2: Variances and covariance for bond and contact num- 
bers at the critical percolation point, divided by N. Nor- 
malizations are chosen such that the curves do not overlap. 
Notice that the covariance is negative for all N and increases 
less fast than N. In contrast, all three variances increase 
asymptotically ~ N, although with vastly different finite size 
corrections. 



Our results for the 2x2 covariance matrix of the bond 
and contact numbers are shown in Fig. 2. More precisely, 
we show there the (co-)variances divided by N, 



C ij = ((ij)-(i)(j))/N. 



(15) 



From this figure we see first that Cbk is negative for all N. 
This is intuitively very plausible: The sum b+k fluctuates 
less than b and k themselves. Moreover, Cbk seems to 
scale as iV a with a w —0.09, i.e. the covariance between 
bonds and contacts increases less fast than TV. Secondly, 
we see that Cbb tends to a constant for N — > oo. A more 
precise analysis shows that the corrections to this are 
oc TV" 1 / 2 , i.e. Cbb ~ (1 -const/iV 1 / 2 ). Next, the variance 
of 36 + k seems to scale with a power of N larger than 
1, but this would contradict Eq. (|13[) . Indeed, a more 
careful analysis shows that the data for Var[36 + k] are 
in perfect agreement with Eq. (|13() , and that the increase 
apparent of Var[36 + k)/N with iV is entirely due to the 
(predicted!) corrections to scaling. Finally, although it 
increases even faster with N, also Var[fc] must ultimately 
scale ~ N, since it is just a linear combination of the 
previous (co-)variances. 

We should stress again that all these results are as ex- 
pected from the scaling theory for percolation (although 
they had not been derived or checked previously, to our 
knowledge |26|). They show that naive power law fits 
without guidance by the scaling theory would lead to 
C k k ~ N with f3 > and thus to erroneous conclusions. 



IV. THE REGION y < 2: COLLAPSING TREES 

Let us first discuss the case y = 0, and treat the gen- 
eral case < y < 2 later. In this case there are only 



0.204 



cvj 0.202 

CD 

m 
in 
d 




FIG. 3: Rfj /N for y = and for five values of r close to 



r = 2, plotted against IniV. 
plot was v = 0.5362. 



The Flory exponent used for this 



the minimal number of bonds, b = N — 1. To find the 
transition point r c = r c (0), we studied the scaling of the 
gyration radius and of the partition sum. We also studied 
the specific heat (i.e., the variance of the number of con- 
tacts), but as in previous cases |2^| this gives much 
less precise estimates. 

In Fig. 3 we plot the re-scaled squared gyration radii, 
R 2 N /N 2l/ , with a suitably chosen value of v, against IniV. 
In Fig. 4 we show effective exponents 9 c g defined by the 
triple ratios |2S| 
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6 In 3 — In 5 

If we assume the scaling ansatz 



(16) 



Z N (y, t) ~ n(y) N N~ e G((r - T c (y))N*) , (17) 

these triple rations should tend to 9 for N — > oo, provided 
t = T c (y). 

From these plots we obtain r c (y = 0) = 2.001 ± 0.001 
(from Fig. 3) and 1.998 ± 0.002 (from Fig. 4). We con- 
jecture that the collapse transition occurs indeed exactly 
at t c = 2. This was already conjectured by previous 
authors [Iflj . although based on much more noisy data. 
Our best unbiased estimates of the critical exponents for 
collapsing trees are 



v = 0.5359 i 0.0003, 9 = 1.842 ± 0.002. 



(18) 



If we accept the conjecture that r c = 2, then these esti- 
mates can be improved to 



v = 0.5362 ±0.0001, 9 = 1.845 ± 0.001. 



(19) 



These values are close to rationals with relatively small 
denominators, which might suggest that v and 9 are in- 
deed simple rational numbers, v — 37/69 = 0.53623... 
and 6 = 59/32 = 1.84375. Together with the estimated 
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FIG. 4: Effective 9— exponents, defined by Eq. H16H . plotted 
against \nN. 
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FIG. 5: Log-log plot of R% /N 2l/ for r = 2 and for five values 
of y (0.0, 0.5, 1.0, 1.5, and 2.0). The Flory exponent used for 
this plot was v — 0.5362 as obtained from Fig. 3. 



TABLE I: Estimates of critical points r c (y), of critical expo- 
nents v and 9, and of inverse critical fugacities fi(y). 
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FIG. 6: Values of ]n[Z N N 9 / y,{y) N ] for r = 2 and for the same 
five values of y as used in Fig. 5. The value of 9 is that for 
y — 0. The critical inverse fugacities [i(y) were chosen such 
that the curves show least bending for large N. 



inverse critical fugacity and with results given below, 
these estimates are collected in Table 1. 

Universality suggests that the same exponents should 
also describe the transition for all y < 2. But direct 
verification is less easy, since there are important correc- 
tions to scaling due to the cross-over from the percolation 
point. Making plots like Figs. 3 and 4 for y > would 
not give very clean results: Searching for cleanest power 
laws would give exponents which depend on y, and would 
give mutually exclusive estimates r c (y) from the gyration 
radius and from the partition sum scaling. 

We therefore adopt a different strategy. Assuming 
that r c (y) — 2 for all y < 2 and that critical expo- 
nents are independent of y, we plotted in Fig. 5 R 2 N /N 211 
against N for different values of y. Similarly, we plotted 
\n[Z^N / fj,(y) N ] in Fig. 6, where fi(y) is carefully chosen 



such as to take into account the dominant exponential 
increase of In Zm with N. In both plots, the curves for 
percolation (i.e. the lower-most curves) become straight 
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[notice that pcrc is usually called r in the 
percolation literature [1|]. The curves for trees (y = 0) 
become horizontal. Finally, the curves for < y < 2 first 
follow the percolation curves and ultimately also become 
horizontal for N — * oo, but very slowly due to the slow 
cross-over. 

Although we cannot claim from Figs. 5 and 6 that this 
scenario is unique, they strongly suggest that there is 
indeed just a slow cross-over from critical percolation to 
collapsing trees, and that the collapse transition occurs 
for all y < 2 exactly at r = 2. 
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FIG. 7: Variances of contact numbers divided by N, plotted 
against r, for y = and for various values of N. 



FIG. 9: Same as Fig. 8, but for number of bonds. 
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FIG. 8: Variances of contact numbers divided by N, plotted 
against ln N. For each y the curve is for our best estimate of 
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FIG. 10: Similar to Figs. 8 and 9, but for the covariances 
between bond and contact numbers. Since these covariances 
are negative for all y and N, we plot actually —G\k- 



Let us finally discuss the fluctuations of the bond and 
contact numbers. The normalized variance of the num- 
ber of contacts, Ckk, of collapsing trees (i.e. at y = 0) 
is shown in Fig. 7 as a function of r for various values 
of N. These data are in good agreement with the less 
precise results of previous simulations 12,14]. They ver- 
ify that r c (0) w 2, but it would obviously be difficult to 
estimate from them r c (0) (as attempted in [T2T[T^) with 
a precision comparable to that obtained from Figs. 3 and 
4. 

Fluctuations for different values of y but precisely at 
t = r c (y) are shown in Figs. 8 to 10. In each figure, the 
normalized (co-)variance is plotted against InTV. For y < 
2 the curves correspond to r = 2, while they represent 
our best estimates of the collapse transition for y > 2. 

As expected, the fluctuations of b are very small and 
tend to zero as y — > 0. They scale ~ N. The bond- 



contact covariances seem to increase slower than N for 
all y > 0, although the decrease of \Cbk\ with increasing 
iV becomes weaker for small y. Of most interest (for the 
present case y < 2) is the scaling of Ckk, since it is the 
contacts which should drive the collapse transition. 

Obviously Ckk increases with N. Naive fits (e.g. 
least square) would give power laws Ckk ~ N 2 ^- 1 with 
4> 0.6 to 0.65. This would agree with previous es- 
timates H, II H3 LI3 i but we should be extremely 
careful with accepting such a fit. The reason is that 
there is no qualitative difference between the cases y < 2 
and y = 2, and for the latter we had seen in Sec. 3 
that Ckk — * const for N — > oo. Thus we propose that 
Ckk — > const, i.e. <j> = 1/2, also for y < 2. If we would 
have <fi > 1/2, then there should be either a clear cross- 
over (which is not seen in Fig. 10), or </> would have to 
be non-universal. 
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V. THE REGION y > 2 

For y > 2 the clusters are richer in bonds, and less rich 
in contacts. The collapse transition was again located 
by requiring both the gyration radius and the partition 
sum to scale. But this time it is much more difficult to 
estimate corrections to scaling. The main reason is that 
our algorithm deteriorates very rapidly when y becomes 
large. While it is still efficient for y < 3.5,jt becomes 
virtually useless for y > 5.5. According to [lfj, the end 
point of the transition line is at y = j/dh = 6.48 ± 0.12. 
At this point, the present algorithm was unable to give 
reasonable statistics of clusters with TV = 200. The algo- 
rithm based on site animals discussed in [23] does a bit 
better and allowed us to obtain a good sample of clusters 
with N — 300. But even for this sample errors are quite 
large, and extrapolation to N — > 00 obviously becomes 
very difficult. We thus do not present any data for the 
Derrida-Herrmann model (r = 0, y = 2/dh), but we just 
state that our simulations were in complete agreement 
with the results of |l6(. 

We therefore must obtain estimates of the critical pa- 
rameters from the region 2 < y < 5, although there could 
again be important cross-over contributions from the per- 
colation point. 

Typical results for y = 3.0 are shown in Figs. 11 and 
12. In Fig. 11 we see that gyration radii suggest the 
transition to be at r c (y — 3.0) = 1.878 ± 0.003, with 
v = 0.5222 ± 0.0005. Effective ^-exponents defined by 
Eq. (|16fl are shown in Fig. 12. They corroborate the 
determination of r c , more precisely they suggest T c (y = 
3.0) = 1.884 ± 0.004. Using as a compromise r c (y = 
3.0) = 1.880 ± 0.003, we obtain as our best estimates 
of the critical exponents v = 0.5220 ± 0.0004 and 9 = 
2.12 ±0.01. 

Estimates of v and 9 compatible with these were ob- 
tained for y = 2.5 and y — 2.25, with r c (y — 2.5) = 
1.970 ± 0.004 and r c (y = 2.25) = 1.994 ± 0.002. But we 
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encountered problems when going to 3.5 < y < 4. Re- 
sults analogous to Figs. 11 and 12, but for y — 3.75, are 
shown in Figs. 13 and 14. This time, there is no value of 
t at which R 2 N shows a pure power law for large N (say, 
N > 100). For small r the curves bend upward (the 
clusters are extended), while for t > 1.630 they bend 
down for very large N, showing that the collapse sets 
first in for very large clusters only. Moreover, taking the 
curve in Fig. 13 for r = 1.625 (which seems to become 
straight for very large A^) to estimate v, we would ob- 
tain v = 0.527, which is definitely much larger than the 
estimates obtained from 2.5 < y < 3.0. 

Similarly, also the partition sum shows late scaling, 
with all curves in Fig. 14 decreasing. Only the curves for 
t > 1.64 shows some tendency to level off for very large 
N. If this is taken as an indication that scaling should set 
in similarly late also for the gyration radii, then indeed 
the lowest curve in Fig. 13 should be the critical one, and 
its slope for N > 1000 should give the correct estimate 
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of v: v ~ 0.520. But all this looks very unconvincing. In 
addition, the value of 9 suggested from Fig. 14 is much 
smaller than that from Fig. 12, any reasonable estimate 
based on Fig. 14 being 9 < 2.06. 

Similar results were obtained for y — 3.5 and 4.0. 
In all these cases the estimated value of 9 is definitely 
smaller than the estimate obtained from y — 2.5 and 
3.0. Our best estimates from the region y > 3.5 are 
9 < 2.0, v < 0.515. Notice that this is not easily un- 
derstood from a cross-over from percolation, since there 
9 = 187/91 = 2.054945. Why should the estimates first 
move away from the percolation value as we move away 
from the percolation point, just to come back to it later? 

Notice that these difficulties, serious as they are for 
the estimation of the critical exponents, have very little 
effect on the estimation of T c (y). In spite of them, our 
estimates (shown in Fig. 1) have errors less than 0.02. 

They have also very little influence on the estimates 
of the bond and contact variances. These depend rather 
weakly on r (thus they would not be very useful for de- 
termining the collapse transition curve, as proposed in 
all previous works). The errors in the data shown in 
Figs. 8 to 10 are comparable to the thickness of the lines. 
As expected, Ckk decreases as y gets larger, while Cbb 
increases. Thus we can safely say that 

Ckk — > const for N — * oo , (20) 

since this was already shown for y = 2. The scaling 
of Cbb is much more subtle. Previous analyses, starting 
with the seminal work of |l6(, obtained Cbb ~ 7V 2< ^ _1 
with 4> ~ 0.6 to 0.65. Just making least square fits to 
the curves in Fig. 9 for the largest values of y would 
give the same estimate. But we think that this is mis- 
leading, as it was also for y < 2. First of all, we would 
expect some sort of cross-over in Fig. 9, which is not seen: 
The transition from y — 2 to y > 2 is perfectly smooth. 
Secondly, we should expect Ckk and Cbb to scale in the 
same way. These two variances control the divergence 



of Zpf(y,T) as we cross the transition curve vertically 
and horizontally. In both cases one crosses the transition 
curve transver sally, and along each transverse line the di- 
vergence should be the same. Therefore we suggest that 
Cbb ~ const, i.e. <j) = 1/2. 



VI. COLLAPSED ANIMALS: ONE OR TWO 
PHASES? 

According to 0, 0, [2^ there are indeed two collapsed 
phases (one rich in contacts, the other rich in bonds). 
The transition line between them bifurcates off from the 
collapse curve not at y = y c = 2, but at y* w 4 (see 
the figure shown in 0]). If this is basically true, but 
with y* w 3.2, this would easily explain our results dis- 
cussed in the previous section: The point (j/*,t*) would 
separate an intermediate part 2 < y < y* of the col- 
lapse curve which is still basically contact driven from 
the bond driven part y > y*. The only unusual feature 
would be that the fixed point (y* , r*) has to be attractive 
when approached from the left along the collapse transi- 
tion curve, while it has to be repulsive on the right. Such 
a behaviour is not typically expected from the RG group, 
although we see no reason why it should be forbidden (it 
requires that the linearized RG flow near the fixed point 
is zero, and the flow is dominated by quadratic terms). 

For a direct test of this scenario, one has to simu- 
late deep in the collapsed region and verify that there 
is a transition between two different collapsed phases. 
To search for such a transition, we have to modify our 
strategy. We cannot use the gyration radius as an in- 
dicator (because R 2 N ~ N in both phases), and we 
cannot use the partition sum either. We checked that 
the free energy of collapsed animals has a surface term, 
lni>jv = Nhxfj, — constiV 1 / 2 + . . ., as one expects for 
any 2-d compact object (a similar term is well known for 
droplets, and was shown to exist also for collapsed un- 
branched polymers |30l l31| ). Even if there is in addition 
a term ~ ln AT in the free energy, there would be just too 
many unknown parameters to pin it down precisely and 
to use this for location the transition. 

Thus we have to look directly to the order parameter, 
which is the difference b — k, and its fluctuations [32|. In 
Fig. 1151 we show the normalized average difference (k — 
b)/N between the numbers of contacts and bonds against 
t, for several values of N. More precisely, we used b — 
N + 1 instead of b in this figure, which is the number of 
bonds exceeding the minimal number. For this figure, y 
was kept fixed at y = 3.75, i.e. we scan Fig. 1 along a 
vertical line which starts at the collapse transition line at 
a r which the previous section suggested to be above r*. 
We see that this difference increases with r, indicating 
that we move indeed from a situation rich in bonds to a 
system rich in contacts. But this increase is so smooth 
that this cannot be used as an argument in favour of a 
thermodynamic phase transition. Instead, Fig. 1151 could 
suggest that there is no such transition. But we should 
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be aware that an analogous search for the 8 collapse in 
ordinary (unbranched) polymers, based on the number 
of contacts, would also fail, because this number varies 
rather smoothly with temperature |33j . Thus we would 
need much larger systems to obtain a firm conclusion 
from a plot like In Fig. I*P31 

Things change if we look at the corresponding vari- 
ances, see Fig. ^| If there is a phase transition, we ex- 
pect a peak which becomes higher as N increases. What 
we see is not quite a sharp peak (it is rather a broad 
bump), but its hight definitely increases with N. Also, 
its maximum shifts with increasing N to higher values of 
r, indicating that it is not related to the collapse tran- 
sition. Notice that the existence of this bump and its 
increase with N is definitely not due to any statistical 
fluctuations. 

It is not clear whether the data shown in Figs. IT51 
and 1161 taken by themselves speak more for or against 
are phase transition between two collapsed phases. They 
certainly do not exclude this possibility. And taken to- 
gether with the anomaly found in the previous section, 
we believe that such a transition is the most natural in- 
terpretation. In this case, a naive (least square, say) fit to 
an ansatz Ck-b,k-b ~ N 2 ^' 1 would suggest <j> ~ 0.54, but 
there is the same kind of curvature in plots of Ck-b.k-b 
versus IniV as in Figs. |H] and [51 and a more careful ex- 
trapolation suggests again that is close to 1/2, more 
precisely (j> = 0.52 ± 0.02. Since the shift of the max- 
imum to larger t with increasing N seems to stop for 
N > 400, we locate the transition near the maximum 
observed for N = 1600, t « 2.1 ± 0.2. Because these 
simulations were very costly in terms of CPU time, we 
have not made similar searches at different values of y, 
and the location of the transition line indicated in Fig. 1 
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should only be taken as a very rough guess. 

We should however point that it seems hard to come by 
with a theoretical argument for such a transition. Usu- 
ally, a phase transition requires some sort of coopera- 
tivity, i.e. some mechanism leading to positive correla- 
tions. In the ferromagnetic Ising model, e.g., this is the 
spin-spin interaction which tends to make neighbouring 
spins parallel. In the antifcrromangetic Ising model on a 
square lattice, the spin-spin interaction favours parallel 
next-nearest neighbours. In the present case we see no 
such interactions. The only source for correlations is the 
fact that the bonds must make the cluster connected, i.e. 
they cannot all be in one part of the cluster, leaving the 
rest of the cluster to contacts. But this should give rise 
to negative correlations, and it is not clear how it can 
lead to a phase transition. 

We should also mention that the individual bond and 
contact variances variances show no bump or peak when 
going into the collapsed region, but fall off continuously 
with r. The bump seen in the variance of k — b is entirely 
due to an increase of the absolute value of the covariance. 
Indeed, as we move off the collapse curve into the col- 
lapsed phase, the scaling of Ckb changes from decreasing 
with N to increasing with it, and continues to increase 
with N even for very large r. 



VII. CONCLUSIONS 

We have presented extensive Monte Carlo simulations 
of collapsing lattice animals and collapsing lattice trees. 
We used a novel algorithm which should be most efficient 
near the (bond) percolation point, and it indeed was. At 
the percolation point it was even more efficient than the 
straightforward Leath algorithm. It was also very effi- 
cient along the entire transition line for y < 2, including 
the point y = which describes the collapse of weakly 
embeddable trees. For the latter we conjectured exact 
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values for all critical exponents. 

Our simulations encountered severe problems for the 
bond-induced collapse at y > 2. For very large y, near 
the point studied first by Derrida and Herrmann, it es- 
sentially breaks down. Thus we were able to verify their 
analysis, but we were not able to improve on it. 

For moderately large y, i.e. for 3.5 < y < 4.0, our 
algorithm is still efficient enough to generate high statis- 
tics samples of rather large clusters (N w 1000). But 
we found there very large corrections to scaling which 
prevent us from extracting precise values of the critical 
exponents, and we can only give upper bounds for v and 
9: v < 0.515, 8 < 2.0. These corrections seem to be 
absent for even smaller y (i.e. 2.5 < y < 3.0), but we 
cannot of course exclude that they would show up if we 
would go to even larger clusters. Assuming this not to 
happen, we can then propose v — 0.5220±0.0004, = 0.5 
and 9 — 2.12 ± 0.01 along the bond-induced collapse 
transition, while v — 37/69 « 0.53623,0 = 0.5 and 
9 = 59/32 = 1.84375 for collapsing trees. Together with 
the exact results v = 48/91 w 0.52747,0 = 1/2, and 
9 = 187/91 « 2.054945 for critical percolation, these val- 
ues show that there are indeed (at least) three different 
universality classes. We cannot rule out that the anoma- 
lies observed at y « 3.75 hint at yet another universality 
class. 

We found somewhat weaker but still statistically sig- 



nificant evidence for a transition between two collapsed 
phases, one contact-rich and the other bond-rich. This 
line branches off from the collapse line at y* = 3.2 ± 0.2, 
which is clearly larger than the value y — 2 for critical 
percolation, but is substantially smaller than the esti- 
mate y* w 4 given in 0, 0, |22| . The cross-over expo- 
nent at this transition is again ss 0.5. Other critical 
exponents were not measured for this transition. 

We should finally comment on our claim that 0=1/2 
for critical percolation (and indeed along the entire col- 
lapse line and for the transition between the two collapsed 
phases). This is based on the definition of in terms of 
bond and/or contact fluctuations. It is also true for site 
percolation, if we look at fluctuations in the number of 
non-wetted surface sites. It should not be confused with 
the fact that the cross-over exponent in a scaling ansatz, 
Eq. 10, usually called a in the percolation literature, is 
smaller than 1/2. From this scaling ansatz it follows sur- 
prisingly that the variances of bond and contact numbers 
must scale ~ N. Our claim that = 1/2 is entirely based 
on this observation, and is made in spite of the fact that 
naive fits would give > 1/2. This should again be a 
warning that power law fits not guided by a solid theory 
can be very misleading. 
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